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The charge and spin Hall conductance of the two-dimensional hole gas within the Luttinger model with and without 
inversion symmetry breaking terms in a perpendicular magnetic field are studied, and two key phenomena are predicted. 
The sign of the spin Hall conductance is modulated periodically by the external magnetic field, which means a possible 
application in the future. Furthermore, a resonant spin Hall conductance in the two-dimensional hole gas with a certain 
hole density at a typical magnetic field is indicated, which implies a likely way to firmly establish the intrinsic spin Hall 
effect. The charge Hall conductance is unaffected by the spin-orbit coupling. 
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INTRODUCTION 

The intrinsic spin Hall effect (SHE) was predicted by 
Murakami et al. 1] in p-doped bulk semiconductors of a 
Luttinger Hamiltonian and by Sinova et al.^ in two di- 
mensional (2D) electron systems with Rashba spin-orbit 
coupling, and the intrinsic SHE has aroused intensive the- 
oretical studies|l|i,[l|l|li,|i[ll[ll[ll[ll[ll[ll[l3 
for its potential applications in the emerging field of 
spintromcs0, 0, 0|. Moreover, experimental obser- 
vation of the SHE has also been recently reported in 
an electron doped sample with the use of Kerr rota- 
tion microscopy [23, |^ and in a two-dimensional hole 
gas (2DHG) by angle- resolved polarization detect ion [2^. 
Kato et 20] suggested that the extrinsic 1^ SHE in n- 
type semiconductor was dominant for that no marked 
crystal direction dependence was observed in the strained 
samples, while the SHE observed in the 2DHG with spin- 
orbit couplin g wa s interpreted as intrinsic by Wunderlich 
et g/.p. la. la l22l|. Further experimental and theoretical 
work is needed to firmly establish the intrinsic VS extrin- 
sic nature of the SHE. Besides, how to develop spintronics 
devices [Hill m is still a very challenging problem. 

Although the SHE in the 2DHG has been intensively 
studied theoretically, most studies have been limited to 
zero external magnetic field 0, H i, 0, H i, ES El ■ The 
spin Hall conductance in the 2DHG arises from the split- 
ting between the light and heavy-hole band and the struc- 
tural inversion asymmetry (SIA)|3,I3| in the 2DHG band. 
In this paper, we study transport properties of the 2DHG 
within the LuttingerP,|23 model with and without inver- 
sion symmetry breaking terms in a perpendicular mag- 
netic field, and two key phenomena are reported. First, 
the sign of the spin Hall conductance {i.e, the direction of 
the spin Hall current) can be modulated periodically by 
the external magnetic field, which means a possible ap- 
plication in the future. Second, a resonant spin Hall con- 
ductance can be observed by adjusting the sample param- 
eters and magnetic field, and the resonance may be used 
to firmly establish the intrinsic SHE in experiments 0|. 
Meanwhile, these two types of spin-orbit coupling have 



no effect on the charge Hall conductance. The rest of the 
paper is organized as follows: The model Hamiltonian 
and the formula for charge and spin Hall conductance 
are introduced in section II. Our numerical result and 
discussion are shown in section III, and the paper is con- 
cluded with a summary in section IV. 



THEORETICAL FRAMEWORK 

The effective Hamiltonian we consider for a 2DHG is a 
sum of both Luttinger and spin-S*^ —3/2 SIA termsP.l23j|: 



2m 

'^ — {Px,Py}{Sx, Sy} + a{S X p)- z 
m 



+ pIs' + {pI)s',) 



(1) 



where S is spin-| operator and m is the bare electron 
mass, and we define {A, B} = \{AB + BA), p'^ = pi + 
Py + iPz)- addition, 71 and 72 are two dimensionless 
parameters modeling the effective mass and spin-orbit 
coupling around the F point. The confinement of the well 
in the z direction quantizes the momentum on this axis 
and it is approximated by the relation (p^) = 0, (p^) ~ 
{Trh/d)'^ for a quantum well with thickness d. This is not 
the most general Hamiltonian when a magnetic field is 
present. The most general one involves both linear and 
cubic spin Zeeman terms, and we neglect these terms 
in Hamiltonian (1), which is a good approximation for 
GaAslHlil. 

We impose a magnetic field B = Bz by choosing 
the Landau gauge A = —yBx, and then p^ = hk^ + 
—y, where — e is the electric charge. The destruction 
operator [23 a = ^ {px + ipy), and the creation op- 
erator flt = , ^ (px — ip,,) are introduced to describe 

Landau levels, where uj = These operators have the 
commutation [a, a^] = 1. In terms of these operators and 
using explicit matrix notation with S = ^ eigenstates in 
the order Sz — +f, +5, —5, — f, Hamiltonian (1) within 
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a magnetic field can be rewritten as 



Hq = huj 



\ 



Hii i^/Zrjo^ — V372a^^ 

^/2i^20? —2ir]a H22 
— \/372a^ —i^fZ-qa 



\ 



-\/372a' 



with 



H22 



(71 + 72) (a^a + ^) + ^(71 - 272), 
(7i-72)(a^a+^) + ^(7i+ 272). (2) 



where the dimensionless parameters 77 = Ci'^\f^^^^ a-nd 

Depending on the confinement scale d the Luttinger 
term is dominant for d not too small, while the SIA term 
becomes dominant for infinitely thin wells, which corre- 
sponds to high junction fields. We first consider the case 



when the confinement scale d is not too small, that is, 
neglect the SIA term: 



2m 

'^ — {Px,Py}{Sx,Sy}. 



The eigenstate of Hq can be taken as 



(3) 



Cnsf24>n-1 
Cnsf3(f>n~2 
\ Cnsf4:(f>n-3 ) 



(4) 



where 0„ is the eigenstate of the -ath Landul level in the 
absence of the Luttinger interaction, and the eigenvalues 
in units of hio when n > 3 can be obtained analytically 
as 



S„+i,/ = (n + i/3)7i + {(71 + 272) + s^-il + [2/(n - /3) + l]7i72 + [(4n - 2{i){n + /) + //? + /^^ + i]7|, (5) 



^+is/ = arctan( 



[71 + /(n - /3)72 + f ] + s^ll + [2/(n - /3) + l]7i72 + [(4n - 2/3)(n + /) + /32 + //? + i]^! ' 



(6) 



r 



The eigenvectors for s = ±1, / = 1 can be expressed as 



^ cos dnsS \ 



sin Qnsf 

V ) 



(7) 



(8) 



while those for s = ±l, / = — 1 can be expressed as 

cos Qrisf 



y sin6'„s/ ) 

where dnsf is defined in Eq.(6), and Qn.\,\ — On, -1,1 + f , 
0n,i,-i = 0n,-i,-i — f • The magnitude of 0nsf is im- 
portant and we will discuss it corresponding to the pa- 
rameters used in this paper bellow, and we can simply 
define that 9n, -i,±i S (0, f ) here. From large n hmit, 
we can deduce that states |n, -|-1,±1) indicate light-hole 
bands (LH^) and |n, — 1,±1) indicates heavy-hole bands 
(HH='= ) y, |25j ■ Besides, the eigenstates and eigenvalues 
for n< 3 are easy to obtain and we do not show them 
explicitly here. The energy levels as functions of the 
fiUing factor v = ^ = are shown in Fig.l. 

We use different color lines to denote different n and 



in order to give a more clear illumination bellow, we 
only plot the energy levels occupied for every v. Solid 
lines indicate Landau levels of mostly holes with spin- 
I (states |n, — 1,1)) and dashed lines indicate Landau 
levels of mostly holes with spin — | (states jn, — 1,— 1)). 
From the bottom up, n increases one by one for solid 
lines (3~8) and dashed lines (O'-^^S) respectively. Param- 
eters used are taken from Ref . . Uh — 2.0 x lO^^/m^, 
71 = 6.92,72 = 2.1, 2^(^)2=40meV, and d=8.3nm. In 
addition, a = 0. Inset: all the parameters used are the 
same except rt/i = 2.73 x lO^^/m^. The dotted line indi- 
cates Landau levels of mostly holes with spin — i. The 
cross indicates the energy crossing which we will discuss 
below. 

Let's begin to study the law in Fig.l in detail (Don't 
include the Fig inset). When i/i — odd integers, the top 
energy levels occupied can be written as 



rpi/i—odd 
occupied 



E. 



-1,-1! 



(9) 



while, when Vi 
occupied are 



even integers, the top energy levels 



TpUi—even 
occupied 



-1,-1,1' 



(10) 



and the total energy levels occupied for every v in- 
clude both states when Vi =odd and Vi =even where 
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FIG. 1: (Color online) Landau levels (units: hu>) as functions 
of filling factors ly. Parameters used are taken from Ref.P, 
and a — 0. Different colors denote different n and only 
energy levels occupied for corresponding i/ are shown. Solid 
lines indicate Landau levels of mostly holes with spin-| and 
dashed lines indicate Landau levels of mostly holes with spin- 
-|. Inset: Uh = 2.73 x lO^'^/m^. The dotted line indicates 
the Landau level of mostly holes with spin- — ^. The cross 
indicates the energy crossing which gives rise to a resonance 
in spin Hall conductance. 



i/i — 1, 2, ...I'. We should add that when n < 3, the defi- 
nition of |n, s, /) is not exact, so we simply take |2, —1, 1) 
as the lowest energy level of n = 2 and |1, 1, — 1) indicates 
the lowest energy level of rt = 1 in the whole paper. It is 
interesting to find that these energy levels arrange regu- 
larly, and the Landau levels of mostly holes with spin-| 
and mostly holes with spin — | appear by turns. This 
law is important and it will determine the behavior of 
the spin Hall conductance which we will discuss below. 

Now we turn to study the system in the presence of an 
electric field applied in —y direction. The full Hamilto- 
nian H is 



(11) 



with 



H = eEy = eE^f^ia^ + a) - (12) 

By treating as a perturbation termls'l up to the first 
order, we calculate the charge and spin Hall conductance 
along X direction, which is the most interesting case. 

The charge current operator for a hole along x direction 
is given by jc = evx with Vx = jfr[x,H], and the spi n 
component current operator can be expressed as[ 



Js =l{S^Vx}. 



(13) 



From Eq.(13) and eigenstates (4) used, we can find that 
the contribution to the charge and spin Hall conduc- 
tance arising from the zeroth order in H is zero for 
< n,s,f\jc,s\n,s,f >= 0. 



We shall only focus on the first order in the perturba- 
tion in H , which is 



{jc,s)nsf 



{n,s,f\jc,s \n',s'J') 
X {n\s'J'\H'\n,sJ)+H.c. 



E 

n's'f 



(14) 



where {jc.s)nsf is the current carried by a hole in the 
state In, s,/ > of , and e„s/ = f>^Ensf- In addition, 
n' — n±l. The average current density of the system is 
given by 



7- '^^ijc,s)nsf fi^nsf), 
nsf 



(15) 



where /(ens/) is the Fermi distribution hmction, and the 
charge or spin hall conductance is 



EL, 



(16) 



We can finally obtain the complete charge and spin Hall 
conductance, and the 2DHG system we consider is con- 
fined in the x — y plane of an area x Ly. 

Our calculation of the charge and spin Hall conduc- 
tance is complete for the linear response theory because 
that the zeroth order mean value of jc,s is zero in each 
state \n,s,f >, meanwhile we estimate the effect of re- 
laxation to the Fermi distribution function fiensf) by 
means of Boltzman equation. The extra term of fiensf) 
from relaxation is proportional to the applied field, thus 
the effect on the conductance from relaxation is second 
order of the applied field E, so we can roughly neglect it 
at linear response region. 



RESULTS AND DISCUSSIONS 

We calculate the charge Hall conductance at T = nu- 
merically firstly. Comparing with the usual quantum Hall 

2 

effect case, i.e, Gc = our results show that the Lut- 
tinger type spin-orbit coupling has no effect on the charge 
Hall conductance, and this result is consistent with the 
quantization of the Hall conductance jl^. 

Now let's turn to discuss the magnetic field dependence 
of the spin Hall conductance. We calculate the spin Hall 
conductance within Hamiltonian (11), and the spin Hall 
conductance as a function of \/B is shown in Fig. 2. In 
order to be convenient to illustrate it. Fig. 2 is divided 
into two regions by a vertical green line at the value of 
magnetic field B = 20T. At the region that B < 20T 
(region 2), the sign of the spin Hall conductance changes 
periodically with the decreasing of magnetic field, and the 
periodicity is ^^^^^ . The spin Hall conductance oscillates 
as a result of the alternative occupation of mostly holes 
with spin- 1 and mostly holes with spin — |. It reaches its 
maxima where i>= odd and minima where i/=even, while 
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FIG. 2: (Color online) Spin hall conductance versus 1/B at 
T — 0, and all the parameters used are the same as those in 
Fig.l. 
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FIG. 3: (Color online) Resonant spin Hall conductance versus 
1/B at T — 0, and all the parameters used are the same 
as those in Fig. 2 except nn ~ 2.73 x Inset: Spin 

Hall conductance in 2DEG with a Rashba spin-orbit coupling 
versus 1/-B at T=0, and data are achieved within the method 
in Ref[ll|. 



the sign of the spin Hall conductance is opposite for odd 
and even v, except when v — 2. 

Let's devote our attention to the region 1, especially 
where ly changes from 3 to 5, and the corresponding B is 
16.5 ~ 27. 6T, which can be achieved in laboratory now. 
When the absolute value of the spin Hall conductance is 
larger than 0.3 (i.e, the spin Hall current is large enough 
to be detected) in the —x direction, the range of magnetic 
field is 20.0 ~ 21. 5T, while in the +x direction, the value 
of magnetic field is in the range of 15.2 ^ 18. 3r. The 
range of magnetic field for a notable spin Hall current is 
wide enough, which makes it easy to modulate the direc- 
tion of the spin Hall current. If the direction of the spin 
Hall current can be detected, we can take the direction 
of the spin Hall current as a new sign, i.e, +x direction 
means "0" and —x direction means "1", and these two 
states maybe the basic for a new logic electronic-device. 
In addition, we have made an scan on hole density, as well 
as magnetic field dependence of sign changes. We find 
that such rich sign changes are robust even in moderate 
magnetic fields, however, a more notable sign changes is 
sensitive to the hole density The sign issue should 
be regarded as a very important aspect of SHE and can 
be used in future experiments (251 i27i| . 

Bernvig et al. 4] have calculated the spin Hall conduc- 
tance within the same Hamiltonian (3) in the absence of 
a magnetic field, and their result shows that the spin hall 
conductance is -^g^ , which is in good agreement with the 
numerical estimate by Wunderlich et a/.l22j. As shown 
in Fig. 2, the jump in the spin Hall conductance tends to 
a stable value of which is very near to the universal 
value. 

In the 2DEG with a Rashba spin-orbit coupling in a 



perpendicular magnetic field, the resonant spin Hall con- 
ductance stems from energy crossing of different Landau 
levels near the Fermi level Although this kind of 
energy crossing can appear when Uh = 2.0 x 10^^/m^ 
of the sample we consider, but energy levels e„=i,i,_i 
and en'=n+i.-i,i are higher than Fermi level as shown 
in Fig.l, so they do not contribute to the spin Hall con- 
ductance, and then there's no resonant spin Hall conduc- 
tance when rih = 2.0 x 10^^/m^ as shown in Fig. 2. If we 
modulate the hole density Uh, and take a value when the 
energy crossing occurs near the Fermi level as shown in 
Fig.l inset {nu — 2.73 x 10^^/m^), there shall be an effec- 
tive energy crossing between en=i,i,-i and e„'=„+i__i^i 
which occurs at Br = 28.52T (the key point denoted by 
a cross in Fig.l inset), we shall see a resonance in the 
spin Hall conductance as shown in Fig. 3. 

In order to make the energy crossing to be occurred 
near the Fermi level, it should request 



3<^<4, 
eBr 



(17) 



that is, when 2.07 x lO^^/m^ < Uh < 2.75 x lO^^^/m^, 
the "effective" energy crossing between mostly holes with 
spin — i and holes with spin-| shall give rise to a reso- 
nant spin Hall conductance in direction as shown in 
Fig. 3. The resonant spin Hall conductance means spin 
accumulation [isl UH. Esj might be observed in experi- 
ments near the edges of a semiconductor channel, and it 
may be a way to distinguish the intrinsic from the ex- 
trinsic SHE©]. 

Even all the results above are obtained numerically, 
we can discuss them in some analytical way. The 
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states near Fermi surface are mostly |n, — 1,±1 >. 
Let's argue the contribution for the spin Hall con- 
ductance when 1/4=0(1(1 firstly, and the key states are 
|n, — 1,— 1 >. We can learn from the properties of 
eigenvectors shown in Eq.(7) and (8) that only states 
|n±l,±l,— 1 > shall contribute to the spin Hall con- 
ductance for {n' ,±1, +1\ H' |ri, — 1,— 1) = 0, furthermore, 
|e„,_i _i - e„±i, -1, -1| < |e„,_i _i - e„±i,i,_i|, so the 
spin Hall conductance arising from states \n,—l,—l > 
shall be dominated by states |n± 1,-1,-1 > 

^ {n, J, \n\s',f) 

, , ^n. — 1, — 1 ^n' s' f 

n's'f 

X {n',s' J' \H'\n, -1,-1) +H.C. 



ddi 



(18) 




0.4 0.6 
1/B(1/T) 



1.0 



and then all the part for Vi—odd can be obtained as 



odd/ 



Yl-K 



9s 



-i), (19) 



odd/ 

9s [n^n) 



where the Landau degeneracy factor muj/{2TiK) is also 
included and 

QB (n,-l,-l|js |n',-l,-l> 

Ech e,i,-l,-l — En', -1,-1 

(n', -1, -1| H' |n, -1, -1) + H.c 
6B (n,-l,l|j, |n',-l,l) 

Ech En, -1,1 — En', -1,1 

{n', -1, l\ H' \n, -1,1) + H.C. (20) 



even 

9s 



(n —)■ n) = 



X 



In deriving Eq.(18), we have used the fact that the tran- 
sitions |i,— 1,— 1 >^ \i — 1,-1,-1 > is canceled by 
\i — 1,-1,-1 >^ |z.— 1,— 1 > where i < n, so only 
the g°g''^{n — > n + 1) is reserved and we simply write 
it as gf'^{n + 1). From Eq.(9), we can obtain that 
ri+l = |- + | + l. In the similar way, all the parts 
for Vi—eyen can be written as 



[Y|^5r"(y)/(e^,-i,+i)].,>4, (21) 

and G^;=2 — ^j- can be calculated immediately. At 



zero temperature, when occupied, any /(e 
The whole spin Hall conductance can be expressed as 

u=odd ~ -r7ry9s \- 

iZTT 



occupied } 



1. 



)+9r{^ + ^)], 



G™„ « Yl^[5r"(^)+5g''(^+3)]. (22) 

The magnitude of is important. Correspond- 

ing to the parameters used, and for the energy levels 
occupied of every ly, we can define that 6'„^_i^i € (0, j) 
and 0„_i_i g (jjf)- From the numerical calculation, 
we can obtain the following linear dependence of filling 
factor that when — odd, 



9rHn) 
n 



1.12n- 2.07,5^"'="(to) w ^ 
i//2 + 5/2,TO= {iy-l)/2 



T.12m-0.97, 

- 1 (23) 



FIG. 4: (Color online) Spin Hall conductance versus 1/B at 
T = 0, Parameters used are the same as those in Fig.l except 
a = 1.0 X lO^m/s. 



and when v = even, 

(n) « 1.12n - 2.04, 5^™"(m) «-1.12n- 0.88, 
n = (i/-l)/2 + 5/2,m = iy/2-l. (24) 



odd 



9s 



From Eqs.(7) and (8), we can also learn that g^'^'^ arises 
from the transition between mostly holes with spin — | 
and i?™*^" arises from the transitions between mostly 
holes with spin-|. The oscillations on the spin Hall 
conductance due to the alternative occupation of mostly 
holes with spin-| and mostly holes with spin — | can be 
explained approximately, and our analytical process also 
indicates that the contribution from the transitions be- 
tween mostly holes with spin-| (or spin — |) should be 
the dominant part of the spin Hall conductance. 

The contribution to the resonant spin Hall conduc- 
tance can be obtained as 



gr = [3V2(7i -t- 72) cos 6*2,-1,1 
•\/2cos02 -1 1 



2\/3- 



72 sm 6'2,-i,iJ 



E2 



El. 



(25) 



-V672 



where ^2 -1 1 = arctan , -, 

' ■ [7l + (f-/3)72]-V[7i + (i-/3)72]^+67|' 

and -©2,-1,1 < which leads to a resonant spin 

Hall conductance in —x direction. A more detailed dis- 
cussion will be given in our further work[25j. 

The relatively large 5 meV measured spitting 0,|22 of 
the HH band implies that the effect of SIA terms is im- 
portant. We have computed the spin conductance within 
the Hamiltonian (2) in the presence of an electric field in 
—y direction, and the result is shown in Fig. 4. 

When ly > 7{i.e, B < 11. 8r), the spin Hall conduc- 
tance oscillates as a result of the alternative occupation of 
mostly holes with spin — |, spin-| and spin — | between 
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some turning points, v =7, 10, 37, 56, 76, 96, while the 
spin Hall conductance oscillates due to the alternative oc- 
cupation of mostly holes with spin-| and spin — | when 
B > 11. 8T. Furthermore, the turning points appear pe- 
riodically when B < 8.3T, and the period approaches a 
stable value 19— V when B < 2.2T. The difference be- 
tween Fig. 2 and Fig. 4 can be found from the effect of SIA 
term. The SIA term arising from the structure inversion 
symmetry rearrange the energy levels when v > 7|2!tI|. 
and this effect is similar to the case in the absence of 
magnetic field |^llo|. 

The charge Hall conductance has also been studied 
within the Hamiltonian (2) in the presence of an electric 
field in —y direction, and our result shows that the SIA 
term has no effect on the charge Hall conductance. Our 
further calculation shows that we can also predict a res- 
onant spin Hall conductance within the Luttinger model 
including the SIA term at a typical magnetic field^^ by 
adjusting the hole density. 

After the completion of this paper, we become aware 
of the independent work 29J of M. Zarea and S. E. Ulloa, 
which studies the system of 2D heavy holes described by 
a fc-cubic model in a magnetic field. In the thin quantum 
well limit, LH bands become energetically irrelevant and 
the HH bands can be effectively described by the 2D HH 
system |^ Ull Isoj where only the lowest HH subband 
is occupied. This condition can be satisfied when the 
quasi 2D system is sufficiently narrow and the density 
and the temperature are not too highj^, 0, Is^l- Al- 
though both the paper by M. Zarea et al. and this paper 
by us compute the spin hall conductance in 2DHG in a 
magnetic field, M. Zarea et al.^^ focuses on the low field 
regime ;30]. In strong magnetic field regime and at a high 
hole density, the resonance reported in this paper can not 
be achieved within the two dimensional HH svstem[25| . 
which is due to the interplay between mostly holes with 
spin — i and holes with spin-|. Moreover, the resonance 
remains within Luttinger model even including Zeeman 
splitting, and the resonance effect stemming from energy 
crossing of different Landau levels near the Fermi level 
due to the competition of Zeeman energy splitting and 
/c-cubic spin-orbit coupling is out of reach for either real 
materials or in experiments |25^. However, when only the 
lowest HH subband is occupied, the 2D HH system can 
reproduce qualitatively the result within Luttinger model 
with Rashba spin-orbit coupling in the presence of a mag- 
netic field j25i|. 

SUMMARY 

In summary, we have studied the charge and spin Hall 
conductance in p-type GaAs quantum well structure de- 
scribed by a Luttinger Hamiltonian with the SIA term in 
a perpendicular magnetic field. The sign of the spin Hall 
conductance (i.e, the direction of the spin Hall current) 



purely caused by the Luttinger type spin-orbit coupling 
can be modulated periodically by the external magnetic 
field, which means a possible application in the future. 
The effective energy crossing between mostly holes with 
spin — i and holes with spin-| at a typical magnetic field 
gives rise to a resonant spin Hall conductance with a cer- 
tain hole density, and the jump in the spin Hall con- 
ductance is very near to the universal value of Our 
calculation shows that the hole density, the well thickness 
and the value of the magnetic field for the resonance are 
accessible in experiments. The resonance may be used 
to distinguish the intrinsic SHE from the extrinsic one. 
Although the SIA term rearranges the energy level of the 
Luttinger hamiltonian, the spin Hall conductance is sim- 
ilar to the case caused by pure Luttinger type spin-orbit 
coupling. The dominate contribution to the spin Hall 
conductance is the transitions arising from |n, — 1,— 1 > 
to |rt±l,— 1,— 1 > (between mostly holes with spin — |) 
and I n, 1,-1 > to |n±l,l,— 1 > (between mostly holes 
with spin-|). Our results also show that both Luttinger 
and spin-S' = 3/2 SIA spin-orbit coupling have no effect 
on the charge Hall conductance. 

We thank Professors Shou-Cheng Zhang, Ruibao Tao, 
Shun-Qing Shen, Wuming Liu, Yu Shi and Dr. B. Andrei 
Bernevig for many helpful discussions. 
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